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Abstract. In the present paper we derive two well-known integrable cases of rigid body 
dynamics (the Lagrange top and the Clebsch system) performing an algebraic contraction 
on the two-body Lax matrices governing the (classical) su(2) Gaudin models. The procedure 
■ preserves the linear r- matrix formulation of the ancestor models. We give the Lax repre- 

sentation of the resulting integrable systems in terms of su(2) Lax matrices with rational 
CN| ■ and elliptic dependencies on the spectral parameter. We finally give some results about the 

^ ' many-body extensions of the constructed systems. 
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^ 1 Introduction 

-i— > 

a 

The Gaudin models were introduced in 1976 by M. Gaudin [5] and attracted considerable in- 
terest among theoretical and mathematical physicists, playing a distinguished role in the realm 
of integrable systems. Their peculiar properties, holding both at the classical and at the quan- 
^ ' turn level, are deeply connected with the long-range nature of the interaction described by its 

commuting Hamiltonians, which in fact yields a typical "mean field" dynamics. 

Indeed the Gaudin models describe completely integrable classical and quantum long-range 
spin chains. The original Gaudin model was formulated as a quantum spin model related to the 
Lie algebra su(2) [5]. Later it was realized that such models can be associated with any semi- 
simple complex Lie algebra q [HI [10] and a solution of the corresponding classical Yang-Baxter 
equation [21 [23]. An important feature of Gaudin models is that they can be formulated in the 
framework of the r-matrix approach. In particular, they admit a linear r-matrix structure that 
characterizes both the classical and the quantum models, and holds whatever be the dependence 
(rational (XXX), trigonometric (XXZ), elliptic (XYZ)) on the spectral parameter. In this con- 
text, it is possible to see Gaudin models as appropriate "semiclassical" limits of the integrable 
Heisenberg magnets [26J, which admit a quadratic r-matrix structure. 

In the 80's, the rational Gaudin model was studied by Sklyanin [21] and Jurco [10] from the 
point of view of the quantum inverse scattering method. Precisely, Sklyanin studied the su(2) 
rational Gaudin models, diagonalizing the commuting Hamiltonians by means of separation of 
variables and stressing the connection between his procedure and the functional Bethe Ansatz. 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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On the other hand, the algebraic structure encoded in the linear r-matrix algebra allowed Jurco 
to use the algebraic Bethe Ansatz to simultaneously diagonalize the set of commuting Hamilto- 
nians in all cases when g is a semi-simple Lie algebra. We have to mention here also the the work 
of Reyman and Semenov-Tian-Shansky [22] . Classical Hamiltonian systems associated with Lax 
matrices of the Gaudin-type were studied by them in the context of a general group-theoretic 
approach. 

Vadim Kuznetsov, to whom this work is dedicated, widely studied Gaudin models, especially 
from the point of view of their separability properties [114 fl2l [T3] and of their integrable dis- 
cretizations through Backlund transformations [HI Q3] • In [13] we collaborated with him showing 
that the Lagrange top can be obtained through an algebraic contraction procedure performed on 
the two-body su(2) rational Gaudin model. Such a derivation of the Lagrange system preserves 
the linear r-matrix algebra of the ancestor model, and it has been used as a tool to construct 
an integrable discretization starting from a known one for the rational su(2) Gaudin model [8]. 

The purpose of the present paper is twofold: on one hand we recall the procedure we used 
in [14] to obtain the Lagrange top from the two-body su(2) rational Gaudin model; on the 
other hand we show how the same technique can be used to derive a special case of the Clebsch 
system (i.e. the motion of a free rigid body in an ideal incompressible fluid) starting from the 
elliptic 5u(2) Gaudin model. In the last Section we show how to construct many-body extensions 
starting from the obtained Lax matrices governing the Lagrange top and the Clebsch system. 



2 A short review of $u(2) Gaudin models 

The aim of this Section is to give a terse survey of the main features of su(2) Gaudin models. 
In particular we shall describe them in terms of their (linear) r-matrix formulation, providing 
their Lax matrices and r-matrices. For further details we remand at the references [5l El [TU[ 

EH [13 Eol E21 EH ESI [26] . 

Let us choose the following basis of the linear space su(2): 

.1/0 — i \ .1/0 -1 \ . 1 / -i 



ai = 2 V "i J ' a2 = 2 I 1 ) ' ^ 
We recall that the correspondence 

t™3 • i i 2 3\ • 1 ( — ia 3 — ia 1 — a 2 \ 

( ' ' a ) — -=- 2 {^ + a 2 ia s ) e«u(2), 

is an isomorphism between (su(2), [•,•]) and the Lie algebra (M 3 , x), where x stands for the 
vector product. This allows us to identify R 3 vectors and su(2) matrices. We supply su(2) with 
the scalar product (•,•) induced from M 3 , namely (a, b) = — 2tr(ab) = 2tr(ba^), Va, b € 
su(2). This scalar product allows us to identify the dual space su*(2) with su(2), so that the 
coadjoint action of the algebra becomes the usual Lie bracket with minus. 

The Lie-Poisson algebra of the iV-body su(2) Gaudin models is given by (minus) ffi^su* (2). 
We will denote by {yf } 3 = i, 1 < z < JV, the set of the (time-dependent) coordinate functions 
relative to the i-th. copy of su(2). Consequently, the Lie-Poisson brackets on 0^5U*(2) read 

{Vi, Vj } = -<Sy e "/3 7 vL 



with 1 < i,j < N. Here e aj a 7 is the skew-symmetric tensor with £123 = 1. The brackets (I2.ip 
are degenerate: they possess the N Casimir functions 

Ci = \{yi,yi), l<i<N, (2.2) 
that provide a trivial dynamics. 
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The su(2) rational, trigonometric and elliptic Gaudin models are governed respectively by 
the following Lax matrices defined on the loop algebra su(2)[ A, A ]: 

N a N 

£ r gW =°«P a + Y, = P + £ x^X^ (2 ' 3) 

i=l 1 i=l 1 

N 1 

4( A ) = £ ■ n n [ ai vi + ff2 y * + cos ( A - Ai ) ^ y < ] ' ( 2 - 4 ) 

sin A — Ai) 
N 1 

£ e( A ) = £^ rr [dnCA-AiJai^ + oij/f + cnCA-Ai)^!/?] , (2.5) 

where the Aj's, with Aj ^ A&, 1 < i,k < N , are complex parameters of the model. We remark 
that in equation (|2.5j) cn(A), dn(A), sn(A) are the elliptic Jacobi functions of modulus k. In 
equation (|2.3p p is a constant vector in M 3 . Its presence is necessary in the rational case in order 
to get a sufficient number of functionally independent integrals of motion. 

It is well-known that the Lax matrices (|2.3p . (|2.4p and (|2.5p describe completely integrable 
systems on the Lie-Poisson manifold associated with ® N su*(2). In particular they admit a linear 
r-matrix formulation, which ensures that all the spectral invariants of Cg(X), Cg(X), Cg(X) form 
a family of involutive functions. Let us give the following result. 



Proposition 1. The Lax matrices C r g(X), Cg(X), Cg(X) given in equations (12. 3p . (12. 4h and (12. 5h 

satisfy the linear r-matrix algebra 

{Cg t,e (X) ® 1, 1 ® £g*' 6 (/u)} + [r rAe (A- / u),4' i ' e (A)0l + l®£^' e (^)] =0, (2.6) 
/or aZZ A, G C, 

r r ,t,e(A) = -/£, e (A) o- Q ® ffa , (2.7) 



and 



f?W = j Va= 1,2,3, 



(/«*).Jfw./fW) i ( s ^j. s ^j.«*w) 

(/;w.j?(A),^w) = f^, 4t 

\sn(A) sn(A) sn(A) 



in equation ()2.6p 1 denotes the 2x2 identity matrix and ® stands for the tensor product in 
C 2 <g>C 2 . 

In the rational case the r-matrix is equivalent to r r (A) = —11/(2 A), where II is the permuta- 
tion operator in C 2 ® C 2 . 

The complete set of integrals of the su(2) rational, trigonometric and elliptic Gaudin models 
can be constructed computing the residues in A = Aj of the characteristic curve det(£g' t ' e (A) - 
fit) = (or equivalents /i 2 = -(1/2) tr[ (£^' e (A)) 2 ]). The following results hold. 

Proposition 2. The hyperelliptic curve det(£g(A) — //I) = 0, A, fi € C, with C r g{X) given in 
equation (|2.3p . provides a set of 2N independent involutive integrals of motion given by 



N i , N N 



J= l ' J i=l i=l 
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Ci = -{yuYi )■ 

The integrals {H^}fL 1 are first integrals of motion and the integrals {Ci}^L 1 are the Casimir 
functions given in equation (12.21 ). 



Proposition 3. The curve det(£g(A) — fi 1) = 0, A, fi G C, with Cg{\) given in equation (|2.4p . 
provides a set of 2N independent involutive integrals of motion given by 

i U ' h * ' 



The integrals {Hj}^L are first integrals of motion and the integrals {C{\f =l are the Casimir 
functions given in equation (j2.2|) . 

Proposition 4. The curve det(£g(A) — /i 1) = 0, A, fj, £ C, with Cg(X) given in equation (|2.5|) . 
provides a set of 2N independent involutive integrals of motion given by 

ge . ^ dn(A t - Aj) yj y] + yj y) + cn(A, - A-,) yj yj " 



3=1 
N 



sn(A; - Aj) 



#o = E [ fi y ) 51 (A * " A i) + f* ? ^ ^( Ai " A i) + ^ f? 93 (Ai - A,-) ] , 



i,j=l 
i+3 



with 



Pioflli^j fooflli^J foi "ll(,AJ 

and 

i?ere a p(\), a, (3 = 0,1, is the theta functior^, and 6 a p = 9 a p(0), 6' a g = (d/d\)\ = Q6 a p(\). 
The integrals {H?}?L Q are first integrals of motion and the integrals {Ci}fL 1 are the Casimir 
functions given in equation (j2.2|) . 

In the rational case it is possible to select a simple and remarkable Hamiltonian. It is given 
by the following linear combination of the integrals of motion {H^f =l given in equation (12.81) : 

N N N 

i=l i,j=l 1 3 i=l 

x We are using the notation adopted in |26| : 

a p{\) = 6 al3 (X,T) = ^exp 7ri (n + ^ 2 t + 2ni (n+ |) (n + | 
a, /3 = 0, 1, where r is a complex number in the upper half plane. 
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where the r/j's with r/j 7^ rjk, 1 < i,k < N , are arbitrary complex numbers. An interesting 
specialization of the Hamiltonian (|2,9|) is obtained considering 7]%= Aj, 1 < % < N: 



(2.10) 



1 JV JV 

«S = « E<y«yi>+EMp,y- 



«J=1 



i=l 



Proposition 5. The equations of motion w.r.t. the Hamiltonian (|2.10p are given by 



N 



k p + E y? ' yi 



1 < i < TV, 



where y» = dyi/dt. Equations (12. lip admit the following Lax representation: 

4(A) = [^(A),^~ ) (A)] = -[£&(A),^ r ' +) (A)], 
Mffli/i i/ie matrix Cg(X) given in equation (|2.3p and 



^'" ) (A) = E 



JV 



i=l 



Aj y» 

A -A,; 



(A) = A P + ^ yi 



t=l 



(2.11) 



(2.12) 



Proof. A direct computation. 



3 Contraction of su(2) Gaudin models: the two-body case 

In the present section we fix N = 2, namely we consider two-body su(2) Gaudin models. 

It is well-known that the Inonii-Wigner contraction of su(2) ©su(2), i.e. a Lie algebra iso- 
morphic to so (4), gives the real Euclidean algebra e(3) [9]. Let us define the isomorphism 
4> e : su*(2) ffisu*(2) -» su*(2) ©su*(2) by the map 

<Ps ■ (yi,Y2) 1 — ► (m,a) = (yi + y 2 ,e(^iyi + ^2)), (3.1) 

where ^1,^2 6 C, 1^1 7^ v 2 and < e < 1 plays the role of a contraction parameter. In the 
limit e — > the Lie-Poisson brackets on su*(2) ©su*(2) are mapped by 4> £ into the Lie-Poisson 
brackets on e*(3) ^ su*(2) © s R 3 : 

{m a , m P} = -e a ^ m 7 , {m a , a 13 } = -e a ^ a 7 , {a a , a?} = 0. (3.2) 

Obviously, the map (f> £ is not an isomorphism after the contraction limit e — ► 0. The Lie-Poisson 
brackets (j3.2[) are degenerate: they possess the two Casimir functions 

ATi = (m,a), K 2 = -{&,&). (3.3) 

A direct calculation shows that if H(yi,y 2 ) and G(yi,y2) are two involutive functions w.r.t. 
the Lie-Poisson brackets on su*(2) ffisu*(2) then, in the contraction limit e — > 0, the functions 
4> £ (H(yi,y 2 )) and </> e (G(yi,y2)) are in involution w.r.t. the Lie-Poisson brackets on c*(3). 

Our aim is now to apply the contraction map <f> e defined in equation (|3.ip to the Lax matrices 
of the two-body su(2) Gaudin models, i.e. the matrices in equations (12. 3p . (|2.4p and (12. 5p with 
N = 2. To do this a second ingredient is needed: as shown in [14 |, I17 |, fTH ], 120] we have to consider 
the pole coalescence Aj = £ fi, i = 1, 2. This fusion procedure can be considered as the analytical 
counterpart of the algebraic contraction given by the map in equation (|3.1|) . 

A straightforward computation leads to the following statement [HI [2D] . 
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Proposition 6. In the limit e — > 0, the isomorphism (|3,ip maps the Lax matrices (|2,3p . ([2 
and (|2.5|) with \ = EVi, i = 1,2, respectively into the Lax matrices 



£ r (A)=P+y + ^, (3.4) 

£*(A) = — — 7TT \o-\ m 1 + (72 m 2 + cos(A) (73 m 3 ] 
sm(A) L J 

+ . o,., [ cos ( A ) fll + ^2 a 2 ) + 0-3 a 3 1 , (3.5) 
sin (A) 

^ 6 ( A ) = — 7TT f dn ( A ) °"1 ml + °"2 m 2 + cn(A) 0-3 m 3 1 
sn(A) L J 

H ^7- ^ I" cn(A) ai a 1 + cn(A) dn(A) g 2 a 2 + dn(A) g 3 a 3 1 . (3.6) 

sn z (A) 



The Lax matrices given in equations (|3.4p . (|3.5p and (|3.6p describe completely integrable 
systems on the Lie-Poisson manifold associated with e*(3). The remarkable feature of the above 
procedure is that the contracted models inherit the linear r-matrix algebra (|2.6p of the ancestor 
system. The following proposition holds \17 \ [T5] . 

Proposition 7. The Lax matrices C r {\), £*(A), C e {\) given in equations (|3.4p . (|3.5p and ([3.6p 

satisfy the linear r-matrix algebra 

{C r,t ' e (X) 81,1® £ r '*' e Gu)} + [r r , t , e (A - /u),/: rAe (A) 01 + 1® £ r ''' e (/i)] = 0, (3.7) 
/or aZZ A,/i £ C, wii/i r r t e (A) giwen m equation (|2.7p . 



3.1 A Lagrange top arising from the rational su(2) Gaudin model 

Recall that the (3-dimensional) Lagrange case of the rigid body motion around a fixed point 
in a homogeneous field is characterized by the following data: the inertia tensor is given by 
diag(l, l,a), a£l, which means that the body is rotationally symmetric with respect to the 
third coordinate axis, and the fixed point lies on the symmetry axis [T| [4l [14} I22j. 

As noticed in p3] the Lagrange top can be obtained from the two-body rational su(2) Gaudin 
model performing the contraction procedure previously described. 

Let us recall the main features of the dynamics of the Lagrange top (in the rest frame). The 
equations of motion are given by: 

rh = [p,a], a=[m,a], (3.8) 

where m € 1R 3 is the vector of kinetic momentum of the body, a S I 3 is the vector pointing from 
the fixed point to the center of mass of the body and p = (0,0, p) is the constant vector along 
the external field. An external observer is mainly interested in the motion of the symmetry axis 
of the top on the surface (a, a) = constant. For an actual integration of this flow in terms of 
elliptic functions see [7]. 

A remarkable feature of the equations of motion (13. 8p is that they do not depend explicitly on 
the anisotropy parameter a of the inertia tensor [3] . Moreover they are Hamiltonian equations 
with respect to the Lie-Poisson brackets of e*(3), see equation (13. 2p . The Hamiltonian function 
that generates the equations of motion (|3.8p is given by 



/[ = i(m,m) + (p,a), (3.9) 

and the complete integrability of the model is ensured by the second integral of motion I\ = 
(p,m). These involutive Hamiltonians can be obtained by computing the spectral invariants 
of the Lax matrix given in equation (|3.4|) . The remaining two spectral invariants are given by 
the Casimir functions of the Lie-Poisson brackets of e*(3), see equation (|3.3p . 
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Proposition 8. The Hamiltonian flow (|3,8p generated by the Hamiltonian (|3.9p admits the 
following Lax representation: 

t r (\) = [C r {X),M {r '-\\)] = -[£ r (\),M (r > +) {\)], 

with the matrix C r (\) given in equation (13. 4p and 

A4 (r ' _) (A) = ^, X (r ' +) (A) = Ap + m. (3.10) 
A 

Proof. A direct verification. ■ 

Remark 1. Using the contraction map (|3,1|) one can obtain equations (13. 8h directly from equa- 
tions OTTT]) (with N = 2): 

m = yi + y 2 = [p,e(^i yi + ^ 2 y 2 )] = [p,a], 

a = e(i/! yi + z; 2 y 2 ) = [yi + yi,e(^i yi + ^ 2 y 2 )] + 0(e 2 ) [m, a]. 

Performing the same procedure on the Hamiltonian Tig = Ai H\ + A 2 i?2 given in equation (|2.10p 
(with N = 2) and on the linear integral Hf + H2 = ( P, yi + y 2 } we recover the integrals of 
motion of the Lagrange top. We have: 

Kg = ^(yi + y2,yi + y 2 ) - Ci - C 2 + (p,e^iyi + e^iy 2 ), 

being C\ = (yi,yi )/2, C 2 = (y 2 ,y 2 )/2 just Casimir functions. Hence, 
7^i=^i<m,m) + (p,a)=Jl. 

Finally, i?[ + ifj = (PiYi + Y2 ) = (p> m ) = ^2- The same procedure allows one to recover 
the auxiliary matrices A4^ r '^(X) given in equation (|3.10p from the matrices Mg (X) given in 
equation (|2.12p . 



3.2 A Clebsch system arising from the elliptic su(2) Gaudin model 

Let us now consider the Lax matrix given in equation (|3.6|) obtained performing the contraction 
procedure on the Lax matrix of the su(2) elliptic Gaudin model with N = 2. 

A direct computation shows that the spectral invariants of C e (X) are given by the following 
quadratic functions: 

If = -(m,m) - -(a,J3ia), 
i| = -(m,Am)--(a,B2a), 
Jfi = (m,a), K 2 = -(a,a), 

where 

B x = diag(0, k 2 ,k 2 - 1), B 2 = diag(0, 0, k 2 - 1), A = diag(l - k 2 , 1,0). 

Obviously, the choice k = in the integrals (|3.1ip and (|3.12p provides the spectral invariants 
of the trigonometric Lax matrix C (A) given in equation (|3.5|) . Thus the system described 
by £*(A) is a subcase of the one described by C e (X). The quadratic functions (|3.1ip and (|3.12p 



(3.11) 
(3.12) 
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are in involution w.r.t. the Lie-Poisson brackets on e*(3) thanks to the r-matrix formulation in 
equation (|3.7|) . 

Let us now recall the main features of the (3-dimensional) Clebsch case of the free rigid body 
motion (in an ideal fluid) [221 127j . This problem is traditionally described by a Hamiltonian 
system on e*(3) with the Hamiltonian function 

ff=i(m,^m)-i(a,Ba), (3.13) 

where (m, a) € e*(3) and the matrices A = diag(«i, a 2 , a 3 ) and B = diag(/3i, 02, /h) are such 
that the following relation holds: 

01-02 j 02-03 | 03 ~ 01 _ q 

a 3 a i a 2 

namely 

&~/?3 /?3~/?l 01-02 

ai = , a 2 = , «3 = , (-3- 14 ) 

72 - 73 73 - 7i 7i - 72 

for some matrix C = diag(7i, 72, 73). 

Taking into account equations (13,lip - (|3.12|) and (I3.14|) we see that C = diag(0, k 2 , k 2 -l) = B\ 
for the Hamiltonian (|3TTT|) and C = diag(l-/c 2 , 1,0) = A for the Hamiltonian ([5712]) . Thus £ e (A) 
can be considered as the Lax matrix of a special case of the Clebsch system described by the 
Hamilton function (|3. 13|) . 

We now derive Lax representations for the Hamiltonian flows corresponding to the Hamil- 
ton functions (|3.11[) — ()3.12|) . They can be written in terms of su(2) matrices with an elliptic 
dependence on the spectral parameter. 

The equations of motion w.r.t. the integrals If and 7| read respectively 

rh=[a,i?ia], a=[m,a], (3.15) 

and 

m = [Am, m] + [a, B2 a] , a = [Am, a] . (3.16) 
A straightforward computation leads to the following result. 

Proposition 9. The Hamiltonian flow (|3,15|) generated by the Hamiltonian (|3.11|) admits the 
Lax representation: 

£ e (A) = [£ e (A),A^(A)], 

with the matrix C e {\) given in equation f)3.6f) and 

M\(\) = — ttt- [ dn(A) ai 01 + a 2 a 2 + cn(A) a 3 a 3 }. 
sn(A) 

The Hamiltonian flow (|3.16p generated by the Hamiltonian (|3.12p admits the Lax represen- 
tation: 

with the matrix £ e (A) given in equation (|3.6f) and 

M%(X) = — T7~T\ [ cn (A) ai mi + cn(A) dn(A) cr 2 ?n 2 + dn(A) cr 3 m 3 ] 
sm(A) 

{ dn(A) 0-1 01 + dn 2 (A) cr 2 a 2 + cn(A) [dn 2 (A) + sn 2 (A)] a 3 a 3 } . 



sn 3 (A) 
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Remark 2. We note that the "traditional" Lax representations for the Hamiltonian flows 
([3.15|) - ([3.16|) are given in terms of Lax matrices depending rationally on the spectral parameter 
[221 127]. However a Lax representation with elliptic dependence on the spectral parameter for the 
Clebsch system is already known [3]. Hence the novelty of our results consists just in establishing 
of the connection between su(2) elliptic Gaudin models and the Clebsch system. 



4 Integrable chains of interacting tops 

As shown in [17\ \19\ [20] one can construct integrable many-body systems starting with the one- 
body Lax matrices given in equations (|3.4|) . (|3.5|) and (|3.6|) . Such systems describe completely 
integrable (long-range) chains of interacting tops on the Lie-Poisson manifold associated with 
© M e*(3), being M the number of tops appearing in the chain. Moreover they admit the same 
linear r-matrix formulation given in equation (|2.6p [17[ [20] . 

Let us denote with (nij, a^) = (mj,mj,mf, a}, a?, a?) € e*(3) the pair of R 3 vectors associated 
with the z-th top of the chain. Thus the Lie-Poisson brackets on © M e*(3) read 

{mf,m,P} = -Sij e Q/ g 7 m], {mf, a 1 ^} = -5 i:j e Q/37 a], {af, o^} = 0, 

with 1 < i,j < M. The above brackets are degenerate: they possess the following 2M Casimir 
functions: 

Cf ) = (m i ,a i > j cf } = ^(a^a*), 1 < * < M. (4.1) 



According to equations f|3.4j) . (|3.5[) and (|3.6|) we can consider the following Lax matrices 
defined on su(2) [ A, A^ 1 ] : 



Al 



A/(A) = P + E^( A -/^)> ( 4 - 2 ) 

i=l 

M 

4, / (A) = ^4(A-/x l ), (4.3) 
i=i 

M 

C% I (X) = ^2ct(X-fi l ), (4.4) 
i=i 

where the fa's with ^ ^ 1 < i,k < M, are complex parameters of the models. The Lax 

matrix C r M (X) describes a system of M interacting Lagrange tops, called Lagrange chain in [TT] . 
while the matrices C l M {X), C e M {X) govern the dynamics of M interacting Clebsch systems. The 
latter models can be called Clebsch chains. 
The following proposition holds |17[ 120] . 



Proposition 10. The Lax matrices C r M (X), C t M (X), C e M {X) given in equations (|4.2h . (14.3 
and (|4.4|) satisfy the linear r-matrix algebra 

{£^ e (A) 01,1® + [r Plt , e (A - /x), £ r A f (A) ® 1 + 1 ® ] = 0, 

for all X,n G C, wii/t r r -t e (A) groen m equation (|2.7|) . 



We now construct the spectral invariants of the Lagrange chain and of the Clebsch chain with 
k = 0. 
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4.1 The Lagrange chain 



The complete set of integrals of the model can be obtained in the usual way. In fact, a straight- 
forward computation leads to the following statement. 

Proposition 11. The hyperelliptic curve det(£^ / (A) — fit) = 0, A,p 6 C, with C r M (X) given in 
equation (|4.2p reads 



t(p,p) + ^E 



i=i 



+ 



+ 



c 



(i) 



+ 



c 



(2) 



X-fj-i (X-fii) 2 (X-fii) 3 (X-fii) 



where 



M 



Rl = (p,xai) + J2 



SI = ( P, a; ) + - < nii, nii } + ^ 



( a j , a^ ) 



A I 



a,, m 



+ 



in - nj (m -Hjf 



The 2M independent integrals {i?[}*£ 1 and {SJ}^ are involutive first integrals of motion and 
the integrals {C^}f£ 1 and }££i are the Casimir functions given in equation (|4.ip . 



M 



Notice that, as in the su(2) rational Gaudin model, there is a linear integral given by Yl R 



i=i 



M 



^(p,m.j). A natural choice for a physical Hamiltonian describing the dynamics of the model 



i=l 



can be constructed considering a linear combination of the Hamiltonians and {Si}f£ il 

similar to the one considered for the rational Gaudin model, see equation (|2.9|) : 



M 



M 



M 



j=l i=l 'J=l 

If M = 1 the Hamiltonian (|4,5p gives the sum of the two integrals of motion of the Lagrange 
top. Our aim is now to find the Hamiltonian flow generated by TL r M and its Lax representation. 

Proposition 12. The equations of motion w.r.t. the Hamiltonian (|4.5|) are given by 



m, 



[P, a; ] + 



M 

MiP + y^m?)"^ 



/^p + x] m J' ai 



(4.6) 



urei/i 1 < z < M. Equations (|4.6|) admit the following Lax representation: 
C r M (X)=\c r M (X),M A r f\x) 



£ r M W,M { li + \x) 



with the matrix C r M {X) given in equation 



and 



M 



M 



(r,-) 
M 



i=l 



Mi m i + 



X&j 



M 



(r,+) 
M 



(A) = Ap + ^nij 



i=l 



Proof. A direct computation. 
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4.2 The Clebsch chain: the case k = 

The complete set of integrals of motion of the Clebsch chain, with k = 0, is given in the following 
statement. 



Proposition 13. The curve det(£* M (A)-/il) = 0, A, /i G C, with £* M (A) given in equation f|4. 3|) 
reads 



1 M 

-fj? = Hi + - [12* cot (A - /ii) + S| cot 2 (A - W ) 
+ Cf } cot 3 (A - /ii) + Cf } cot 4 (A - Mi)] , 



where 



M M / M \ 2 

^ = i E [ k 1 ) 2 + ) 2 ] - 1 E ™? + M E *n 

i=l ij=l \i=l / 

i+3 

M 1 

+ E sin(/it _^.) [ a i m j + « 2 m l + a ? m ) cos (^ - to) ] 

2,7 = 1 



1 ^ ^ ^Q^if-^i f^i) rll 2 2 / \ i 

+ 2 E sinfc.-^.) L a * G i + a * a J + a * a J COS ^ ~ W\ ' 

i,?' = l 



1 M 



7=1 



+ E ^7 — " T [ m i m ) + m i m i + m ? m j cos(^i — Mi) ] 

sm(/ij — MjJ 



r- 

3+1 
M 



+ ^ -r~, — \ m, a,- + m,- a, + mf af cosf/Ui - Mf) 

^— ' sm(uj — M?) 
i=l 



7= 
3+1 

11 2 2 3 3 

m^- Oj — m^- Oj — c 
M 



(lH-f*j)] 



2 E sinS^.^ ) [ a * a ) + °Z a2 J + ^ a i cos(/li " ^ ] 



7=1 



i 1 1 M 

= cf + 1 [K 1 ) 2 + (m 2 ) 2 + (mf) 2 ] + i(«3) 2 + I J2 «M 



i,j = l 
3+i 



+ E ^(/J- [ a * mjl + a * 2 m2j + ° 3 mj3 cos( ^ " h) ] 



7 = 
7'A 



Ecotf^i^ A^fJ rii 22 / \i 

sin^-^.) I a * a i + a * a 7 + a i a 3 cos ^ - Mj) J ■ 

7=1 

7V* 
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The integrals Hq, {i?-}^, {Sj}f£ 1 are involutive first integrals of motion (only 2M of them 

are independent). The integrals {C^ and {Cf' }f£ 1 are the Casimir functions given in 

equation (|4.ip . 

5 Concluding remarks and open problems 

In the present paper we have proposed an algebraic technique which enabled us to derive two (3- 
dimensional) integrable cases of rigid body dynamics (the Lagrange top and the Clebsch system) 
from two-body su(2) Gaudin models. We remark that the explicit construction of the Lagrange 
top starting from the su(2) rational two-body Gaudin system has been presented for the first 
time in [TJ]. To the best of our knowledge the derivation of the Clebsch system defined by the 
involutive Hamiltonians (|3.1ip - (|3.12p starting from the su(2) elliptic two-body Gaudin system 
is new, although the novelty is essentially in establishing of the connection between these two 
integrable systems. 

Let us stress that the construction outlined here is just a top of an iceberg. In [T71 Q2J HU [20] 
we presented a general and systematic reduction, based on generalized Inonii-Wigner contrac- 
tions, of classical Gaudin models associated with a simple Lie algebra 0. Suitable algebraic 
and pole coalescence procedures performed on the iV-pole Gaudin Lax matrices, enabled us to 
construct one-body and many-body hierarchies of integrable models sharing the same (linear) 
r-matrix structure of the ancestor models. This technique can be applied to any simple Lie 
algebra and whatever be the dependence (rational, trigonometric, elliptic) on the spectral pa- 
rameter. Fixing g = 5u(2), we constructed the so called su(2) hierarchies \18\ I20|. In particular 
the Lagrange top corresponds to the first element (N = 2) of the su(2) rational hierarchy, and 
the Clebsch system is the first element of the su(2) elliptic hierarchy. 

We studied also the problem of discretizing the Hamiltonian flows of the su(2) rational Gaudin 
model. One of the authors (O.R.), together with Vadim Kuznetsov and Andy Hone, constructed 
in [8] one-point (complex) and two-point (real) Backlund transformations (BTs) for this model. 
Later on, in p3], again in collaboration with Vadim, we studied the problem of discretizing the 
dynamics of the Lagrange top using the BTs approach [151 CEE] • 

In [201 121j . using a different approach, we have obtained a new integrable discretization for 
the Hamiltonian flow given in equation (|2,lip . It is expressed in terms of an explicit Poisson 
map and a suitable contraction performed on it enables us to construct discrete-time versions of 
the whole su(2) rational hierarchy. Our results include, as a special case (N = 2), the discrete- 
time version of the Lagrange top found by Yu.B. Suris and A.I. Bobenko in [4]. Moreover, the 
same procedure enabled us to find an integrable discretization of the Hamiltonian flow (|4.6p . 
describing a discrete-time version of the Lagrange chain. 

A natural extension of our discretizations could be the construction of a suitable approach 
for models with a trigonometric or elliptic dependence on the spectral parameter instead of 
a rational one. To the best of our knowledge there are no results in this direction in literature. 
We remark here that integrable discretizations for the flows (|3.15[) — f)3.16[) have been found by 
Yu.B. Suris, see [27l EHl [29] , by using rational Lax matrices. 
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